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Consider a modified, stable, two node Jackson network where 
server 2 helps server 1 when server 2 is idle. The probability of a 
large deviation of the number of customers at node one can be cal¬ 
culated using the flat boundary theory of Schwartz and Weiss [Large 
Deviations Performance Analysis (1994), Chapman and Hall, New 
York], Surprisingly, however, these calculations show that the pro¬ 
portion of time spent on the boundary, where server 2 is idle, may 
be zero. This is in sharp contrast to the unmodified Jackson network 
which spends a nonzero proportion of time on this boundary. 


1. Introduction. In this paper we derive the rough (logarithmic) asymp¬ 
totics for the steady state probability 7r of a particular two node queueing 
network as the queue at server 1 gets large. The analyzed queueing network 
is a variation of a two node Jackson queueing network in which server 2 when 
idle can assist server 1. Allowing one of the servers to help can completely 
change the behavior of the network. This network clearly exhibits a large 
deviation phenomenon, which we call a bridge. For certain parameters, as 
the queue length at node 1 grows, the queue length at node 2 stays small, 
but generally positive so that server 2 is prevented from helping server 1. 
Instead of jittering along the x-axis, the process skims above the x-axis and 
only rarely touches the axis. 

This bridge phenomenon seems to have been somewhat overlooked. In 
particular, the theory in [8, 11] for analyzing exact asymptotics does not 
apply. In a companion paper [9], we extend the theory and develop an ap¬ 
proach to obtaining the exact asymptotics of networks exhibiting the bridge 
phenomenon. 
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The bridge phenomenon in the modified Jackson network is not an isolated 
case—the bridge phenomenon is ubiquitous. Since becoming aware of it, we 
are encountering it frequently in a variety of contexts. In Section 5 we use 
this theory to revisit the bathroom problem discussed by Shwartz and Weiss 

[15]- 

In Section 2 we describe the Jackson network and the modified network. 

We then discuss the possible large deviation paths for overloading node 1. 
Section 3 determines the stability conditions of the modified Jackson net¬ 
work. Section 4 analyzes the rough asymptotics of the modified Jackson net¬ 
work; the analysis of the exact asymptotics appears in [9]. Section 5 briefly 
describes another model where the bridge phenomenon occurs. 

2. Notation and main results. Consider a Jackson (1957) network with 
two nodes. The arrival rate of exogenous customers at nodes 1 and 2 form 
Poisson processes with rates Ai and A 2 , respectively. The service times are in¬ 
dependent, exponentially distributed random variables with mean 1/pi and 1 / ^2 , 
respectively. Each customer’s route through the network forms a Markov 
chain. A customer completing service at node 1 is routed to node 2 with 
probability r\^ or leaves the system with probability ryo := 1 — r\p,. Rout¬ 
ing from node 2 is defined analogously. So without loss of generality, we are 
assuming rip = r' 2,2 = 0. The routing process, service processes and arrival 
processes are independent. 

To ensure that the network is open, we assume that n 2 ^ 2,1 < 1- Since the 
network is open, the traffic equations 

(2.1) A* — A* T \ 3 -ir 3 -ij for i — 1,2, 

have a unique solution (Ai, A 2 ) = ((Ai + A 2 r 2 ,i)/(1 - ri^r^i), (A 2 + Airi j2 )/(1 — 
n, 2 f 2 ,i))- To eliminate degenerate situations, we assume that Ai > 0 and 
A 2 > 0. 

The joint queue length process of this Jackson network forms a Markov 
process with state space S = {0,1,.. .} 2 . Define pt = Xi/p-i, for i = 1,2. From 
Jackson (1957), it follows that the stationary distribution for the joint queue 
length process being in the state (x,y) G S' is (1 — pi)pf (1 — P 2 )p\i provided 
that the stability conditions p\ < 1 and P 2 < 1 hold. 

The network that we analyze is a small change from the above network. 
Suppose that server 2 has been cross-trained and helps server 1 whenever 
queue 2 is empty. Let p\ > p,\ be the combined service effort of the two 
servers at node 1 when server 2 is empty. The transition rates for the 
joint queue length process of the modified network are shown in Figure 1. 

By comparing jump rates, the total number of customers in this modified 
network is stochastically smaller than the total number in the associated 
Jackson network. Hence, the modified network will be stable if the asso¬ 
ciated Jackson network is. However, cross-training server 2 may allow the 
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modified network to be stable even if p\ > 1. In particular, if /02 < 1 and 
p\ > (Ai — PXP 2 )/ (1 — P 2 ), then the modified network is stable; see Section 3 
for the argument. 

We are interested in the rare event of a large deviation in the number of 
customers at node 1; that is, more than £ customers at node 1 where £ is 
large. The steady state probability of this rare event is proportional to the 
number of visits to Fg = {(x, y) : x > £, y > 0} between returns to the origin. 
For the Jackson network, we can determine the most likely path from the 
origin to Fi by looking at the reversed process starting from steady state in 
F( and look at sample paths that leave F( on the first step and never return. 
The reversed process has external arrivals entering node i with rate A^r^o- 
The service rates at the nodes are unchanged, but the routing probabilities 
for the reversed process are rT = Xtrij/Xj for i and j in {1,2}. Start the 
reversed process in state (£,y), where £ is large but y is small. As long as 
there are customers at node 1, customers leave node 1 at rate p\. Thus, 
customers enter node 2 at rate A 2 ?’ 2 ) o + hi r i 2 - If this rate is less than the 
maximum rate at which customers can leave node 2, that is, 

(2.2) A 2 r 2i0 + Pir\p < P 2 or, equivalently, p ^ 1 > r 2 , 0 + rpAP^ 1 , 

then the number of customers at node 2 remains small, and the reversed 
process starting from F( bounces along the x-axis to (0,0). If the inequality 
is reversed, 

P 2 1 < r 2fi+r2,iPi 1 -, 

then the process starting from (£, 0) leaves the x-axis and heads roughly 
northwest (with an easily determined slope) as the customers in node 2 
grow until node 1 empties, that is, hits the y-axis. 
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From there, the process bounces along the y-axis south to the origin 
because customers in the time reversed network enter node 1 at rate Airgo + 
P 2 r 2 .i- This rate is less than the maximum rate at which customers can leave 
node 1, that is, 

(2.3) Airi j0 + ^ 2 ^ 2,1 < Mi or, equivalently, p j -1 > n j0 + r^p^ 1 , 

given (2.2) fails (otherwise just add the two inequalities together and derive 
a contradiction). 

If the inequality in (2.2) is changed to an equality, then the number in 
node 2 in the reversed process behaves like a simple, symmetric random 
walk, which would hit the y-axis at a height proportional to \fl. Thus, in 
the Jackson network, there are three possibilities for the most likely approach 
from the origin to F), though the approach corresponding to equality in (2.2) 
occurs only for a set of parameters with Lebesgue measure zero. 

Now consider the modified network with p\ = p\. Of course, this is identi¬ 
cal to the Jackson network. Suppose the most likely path for a large deviation 
at node 1 of this Jackson network bounces along the x-axis. Let p\ increase. 
As pi increases, the approach going out the x-axis becomes more difficult 
and may eventually become more difficult than some other approach. In 
addition to the obvious possibility of going up the y-axis, it turns out that 
there is a third possibility hinted at by the case of equality in (2.2): the pro¬ 
cess travels along the x-axis, but instead of jittering along the x-axis, the 
process skims above and only rarely touches the x-axis. This third approach 
we refer to as a bridge path , which is slightly optimistic since we hope to 
prove properties in a later paper that would justify the word “bridge.” 

Section 4 contains our preliminary investigation of this modified Jackson 
network by looking at the behavior under the fluid scaling. For the fluid 
scaling, speed up the transition rates by l and measure customers in units 
of 1 /t, which results in a functional s.l.l.n. In particular, we apply the flat 
boundary theory of Schwartz and Weiss [15] to obtain rough asymptotics, 
as well as the fluid scaled large deviation path. In the fluid scaling, both the 
bridge path and the path that jitters along the x-axis collapse to a constant 
speed line along the x-axis, which suggests that the flat boundary theory 
might not be able to distinguish between the two. However, the calculations 
also give the proportion of time spent on the boundary. In some cases, the 
proportion of time spent on the boundary is zero, proving the existence of 
this third possible approach to Fg. In fact, we define a bridge path to be 
such a large deviation path which follows a line, for example, an axis though 
the proportion of time the process spends on the line is zero. Although the 
term is defined with respect to the fluid scaling, the basis for the term is the 
conjectured behavior of the unsealed process. Even though we suspect that 
it has a bridge shape, there are other possibilities. For example, the most 
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likely path when equality holds in (2.2) jitters up the y-axis proportional to 
\fl before drifting to Fg also spends zero time on the boundary and collapses 
to the x-axis under the fluid scaling. We intend to sort out these questions 
in a future paper. 

A jitter path follows a line, for example, one of the axes, while spending 
a nonzero proportion of time on the line. Though jitter path is defined with 
respect to the fluid scaling, the term reflects the behavior of the unsealed 
process, which jitters along the line as it travels to F(. 

We will use the phrase “with large deviation rate 0” to mean that 

lim \\ogP{W G F(\W{fS) = (0,0)) = —0, 

l —>oo t 

where F(£) is the set of cadlag paths in S starting at the origin and associ¬ 
ated with a large deviation of W to F( before returning to the origin; that 
is, to describe the rough asymptotics. Note F(t) is a set of paths hitting the 
set F(_. 

Basically, we will show that the rough asymptotics can be determined 
from three points; see Figure 2. The coordinates of the easternmost point 
of the egg-shaped curve is labelled 9 b . If the curve M~ = 0 intersects the 
egg M + = 0 between 6 b and the 0i-axis, then the intersection is labelled 0 J ; 
otherwise, 0 J = 6 b . If the horizontal line at height log(p 2 1 ) intersects the egg 
between 9 b and the 0i-axis and (2.3) holds, then the intersection is labelled 
0 C ; otherwise, 9 C = 9 b . 

The first coordinate of 9 b , 0 J and 9 C gives the large deviation rate of the 
best bridge path, jitter path and cascade paths, respectively. The minimum 
of the three first coordinates is the rate associated with a large deviation at 
node 1. Theorem 4 summarizes these results. 

3. A bound and stability. We will need the following bound in [9]. Since 
the stability argument and the derivation of the bound use the same cou¬ 
pling, we have included both in this section. 

Lemma 1. For the stable, modified network, 

(3.1) T,*{0,j)<cfr 

j>y 


Proof. The proof is divided into two parts. First, we use a coupling 
argument to consider the case < fi\. We associate a region R(x,y) with 
each point (x,y) E S. The argument shows that there is a coupling of the 
queue length processes of the Jackson and modified networks so that if the 
modified network is in state (x,y), then the state of the Jackson network is 
in R(x, y). It immediately follows that the stationary probability of the mod¬ 
ified network being in state (x, y) is bounded by the stationary probability 
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that the Jackson network is in R(x,y). Similarly, it follows that the station¬ 
ary probability that the modified network is in B C S is less than or equal to 
the stationary probability that the Jackson network is in \Jr x y \ eB R(x,y). 

Now we describe the coupling. We can consider them to be a pair of 
discrete time Markov chains, W[n] for the Jackson network and Y[n] for 
the modified network, subordinated to a common Poisson process with rate 
Ai + A 2 + Hi + / 12 , which we assume without loss of generality to be one. 
Basically, W[n] and Y[n] “attempt” to move in the same direction. More 
precisely, generate an i.i.d. sequence of random variables (directions) taking 
values {E, N, W, NW, S, SE, W*, NW*} with probabilities {Ai, A 2 ,/ri(l — 
^ 1 , 2 ),Mir- 1 , 2 , M 1 — ^ 2 , 1 ),A* 2 ^ 2 ,i, -ri >2 ), (nl~y,i)r lj2 }, respectively. 



Fig. 2. M+ = 0 is the egg-shaped curve; the other curve is M = 0. 
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For the nth step, both W[n] and Y[n] “attempt” to take a single step in the 
direction given by the nth random direction. The directions W* and NW* 
indicate west and northwest, but only for the modified network when it is 
on the x-axis; otherwise, W* and NW* indicate that the process stays put. 
By “attempt,” we mean that the process moves to the neighboring state in 
that direction unless it would result in the process leaving the state space; 
for example, if either process were on the y- axis and the attempted step 
were to be in the direction W. Note that if both processes move, they move 
in the same direction. However, one process may move while the other stays 
put. 

For the regions, define 

fo 9 \ R{x,y) = {( i,j ) <E S\i > x,j > (y - 1)+} 

1 j U{(i, 3 )€S\(i,j) = (x-l,y)}. 

Note that for (x,y) > (1,1), both (x — 1 ,y) and (x,y — 1) are in R(x,y), but 
(x — 1, y +1) and (x — 1, y — 1) are not. We claim that if W[n] G R(Y[n]), then 
W[n + 1] G R(Y[n+ 1]). This is clear when both processes are in the interior 
since they both move in the same direction. Consider the case Y[n] = (1 ,y) 
with y > 1 and W[n] = (0, y). Both processes move in the same direction 
unless the direction is W or NW', in all cases, W[n + 1] will still be in 
R(Y[n + 1]). Now consider the case when Y[n] = (1,0) and W[n] = (1,0). 
Consider the movements NW*, then W and then SE. This makes W[n + 
3] = (0,0) and Y[n + 3] = 1 and W[n + k] will still be in R(Y[n + k]) for 
k = 1,2,3. This trajectory explains why R{x,y) is defined as it is. We leave 
it to the reader to finish checking the claim. 

Now start W[0] and y[0] off in the same state with distribution n. Notice 
that 

E *{i,j) = P{Y[n]€R(0,y + l)) 

i>0,j>y 

= P(W[n] G R(Y[n]),Y[n] G R( 0, y + 1)) 

< P(W[n\ € R(0,y)) < p ( 2 y ~ 1)+ 

for sufficiently large n because W[n] converge to steady state. This completes 
the argument under the condition Ai < y,\. [R(x,y) is the smallest set that 
will work for this coupling. This can be seen by starting both chains at the 
origin and arguing that there is a sequence of random directions such that 
(Y[n], VF[n]) = ( x,y,i,j ) for every ( x,y ) G S and (i,j) G R{x,y). Most points 
can be reached by having both processes move sufficiently far east, then the 
modified network moves west back to the origin, then the Jackson network 
moves sufficiently far northwest, and finally both processes move north and 
east sufficiently far. The few remaining points can be reached by getting the 
Jackson network to the origin while the modified network is at (0,1).] 
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Now consider the case when Ai > (i\. The number of customers at node 2 
is stochastically smaller than a birth-death process on the nonnegative in¬ 
tegers with birth rate A 2 + hi ? ’i ,2 and death rate P 2 on the positive integers, 
and birth rate A 2 + //*ri ,2 in state 0. Thus, the probability of more than y 
customers at node 2 is smaller than c[(A 2 + ^ 1 ^ 1 , 2 )/ p 2 \ y for a suitably chosen 
constant c. Using (2.1) and that Ai > yui, it follows that [(A 2 + Pirip) / p 2 \ < 
P2- □ 

3.1. Stability of the modified network. 

Proposition 1. The joint queue length process of the modified network 
is positive recurrent if X 2 < P 2 and \\ < P 2 P 1 + (1 — P 2 )p\- If either inequality 
is reversed, then the process is transient. 

Proof. The result follows from comparing the modified and Jackson 
networks in the case p\ < 1; for the coupling, see the first part of the proof 
of Lemma 1. Now consider the case where Ai > p\. The following definitions 
of T n and X n are only used in this proof. Let T n denote the time that the 
nth busy period starts at node 2, and let X n be the number of customers in 
queue 1 just prior to the start of the nth busy period at node 2. The process 
Xq,X\, ... is a Markov chain. To prove the result, it suffices to show that 
E[T n+ i — T n \X n ] < 00 and that Xo,X \,... is positive recurrent. 

The random variable T n +i — T n represents the nth busy cycle at node 2, 
which is the sum of the nth busy period and the nth idle period at node 2. 
Since A 2 > 0, the expected length of the nth idle period is finite. The length 
of the nth busy period at node 2 is stochastically increasing in X n \ hence, 
it is stochastically smaller than the busy period in an M/M/1 queue with 
service rate P 2 and arrival rate A 2 + AH r i, 2 - Fortunately, A 2 + hi r i ,2 < P 2 
since A 2 + /DM ,2 < A 2 when Ai > p\. The expected length of this bounding 
busy period is l/(// 2 - (A 2 + Mi r i, 2 ))- 

To show that Xq. X \,... is positive recurrent, we will show that for all X n 
greater than some large constant M, E[X n+ i — X n \X n ] < — e < 0, which is 
a Foster-Lyapunov type condition guaranteeing stability. The constant M 
can be chosen large enough so that the busy period at node 2 is arbitrarily 
close to the busy period of an M/M/1 queue with service rate p .2 and arrival 
rate A 2 + /DM ,2 and that the departure process from node 1 during this busy 
period is arbitrarily close to a Poisson process with rate The change in 
the queue length at node 1, X n+ \ — X n , can be decomposed into the sum 
of the change during the busy period at 2, and the change during the idle 
period at 2. Hence, for X n > M, E[X n+ i — X n \X n \ is arbitrarily close to 

(Ai + ^ 2 i’2,i ~ Pi) (Ai - p\) 

{P 2 ~ (A2 + /DM,2)) (A2 +/LjVi^) ’ 


(3.3) 
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which is strictly less than zero if p*\ > (Ai — p\p 2 )/(l — P 2 ) and is equivalent 
to p 2 pi + (1 - P 2 )P*\ > Ai. 

To show transience, first assume that A 2 < P 2 , but Ai > P 2 P 1 + (1 — P 2 )p\- 
From the argument above, for X n > M sufficiently large, E[X n _|_i — X n \X n \ > 
e > 0; hence, for M sufficiently large, X n behaves arbitrarily closely to a 
random walk with strictly positive drift and is transient. 

Now consider the case when A 2 > P 2 - Note that A* represents the long run 
average arrival and departure rate from node i assuming that the network, 
either Jackson or modified, is recurrent. Assume X n is recurrent; otherwise, 
we would be done. Since X n is recurrent, the departure rate from node 1 
equals the arrival rate. If this rate is less than Ai, then the queue length at 
node 2 must be diverging. However, if the departure rate from node 1 is Ai, 
then node 2 is also diverging since the arrival rate would be A 2 - In either 
case, the network is transient. □ 

4. Flat boundary approach. We are interested in describing how a mod¬ 
ified Jackson network overloads. In particular, we will be interested in the 
rare event when the system starts out empty and node 1 reaches a level i 
before the system empties again. 

If we speed up the jump rates by a factor £, but reduce the jumps by a 
factor \/l, then we get the scaled process W£. From the theory in [15], 

lim \ log P(W £ F(^)|VF(0) = (0,0)) 

£—*00 z 

= lim )logP(W e G F(l)\W e (0) = (0,0)) 

1—>00 £ 

(4.1) = - inf [I(p)], 

peF 

where I(p) = J 0 T A(^|(s),p(s)) ds is the rate function associated with a path 
p in the set of absolutely continuous paths F starting from (0,0) which hits 
the set {(l,y) :y > 0} at some time T before returning to (0,0). Note that 
after time T, the path follows the natural drift path back to the origin, 
and A(^(s),p(s)) = 0 along the drift path. We will use better, cheaper and 
smaller action, synonymously. 

Since W has constant (but different) jump rates on and off the ffat bound¬ 
ary, the local rate function A (v,w) = A + (F) if w is in the interior; that is, 
if w = (x,y) with y > 0 and A (v,w) = A - (i;) if w = (x,0). It follows from 
the calculus of variations that the cheapest path in F is a sequence of line 
segments of constant speed which changes direction only on the y-axis. [If 
a path changes direction on the x-axis, either leaving the interior to travel 
along the x-axis or by leaving the x-axis at some point other than the ori¬ 
gin to travel through the interior, there is a cheaper path from the origin 
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to {(l,y) :y > 0}; the cheaper path might hit {(l,y) :y > 0} at a different 
point like (1,0).] Consequently, the cheapest path in F must lie in one of 
the following sets Fi,F x ,F c , where: 

Fi is the set of all constant speed paths with positive slope across the inte¬ 
rior until hitting (1 ,y) with y > 0, 

F x is a constant speed path jittering or forming a bridge along the x-axis 
until hitting (1,0), and 

F c is a constant speed path jittering up the y-axis, then changes direction 
and heads for the point (1,0) at constant velocity. Thus, the customers 
first build up in node 2, and then cascade into 1. 

In each of these three cases, we will be able to reduce the problem to a differ¬ 
entiable, constrained nonlinear optimization problem. We use the Karush- 
Kuhn-Tucker conditions, which are given in a variety of texts including 
[3, 13], to determine the minimal action in each of the three cases. These 
results are then combined to determine (4.1). 

4.1. A bridge is better than any path through the interior. In this section 
we consider the interior paths T) C F. That is, we consider paths that ini¬ 
tially have the formp(s) = (v\s,V 2 s) until hitting (l,u 2 /ui) at time T = l/v\, 
where v\ > 0 and V 2 > 0. After time T, the path follows the natural drift 
path until reaching the origin. We will find the inf pG [I(p)] and show that 
this inf is not attained in i 7 ). However, there will be a bridge path which 
attains the inf. Hence, there will always be a bridge that is better than every 
interior path. 

Define the log moment generating function of the compound Poisson pro¬ 
cess associated with jumps in the interior 

M + (6 1: e 2 ) = Ai(e ei - 1) + MiD,o(e“ 01 - 1) + Piri, 2 (e 02 ~ 01 - 1) 

+ A 2 (e 02 - 1) + /r 2 r- 2 ,o(e _02 - 1) + P 2 r 2 ,i(e dl ~ 02 - 1). 

The Hessian of M + is positive definite; hence, M + is strictly convex. We also 
know that M + (0,0) = 0. Now, we argue that there exists a point ( 61 , 62 ) > 
(0,0) with M + ( 61 , 62 ) < 0. To see this, first consider the case p\ < 1. In 
this case, VM + (0,0) • (1,1) < 0. Hence, (1,1) is a decreasing direction. Now 
assume p\ > 1. In this case, (0,1) and (e, 1) are decreasing directions for a 
suitably small e > 0. 

Next, the local rate function in the interior in direction v = (ui,u 2 ) is 
A + (v) = sup(0iui + 62 V 2 - M + (0i,6> 2 )), 

61,62 

which is clearly convex. For a path p € F) with velocity v, A(^(s),p(s)) = 
A + (u) for s > 0. Hence, I(p) = A + (v)/v\ is a good rate function; see Theo¬ 
rem 5.1 of [15] . The remainder of this section is devoted to finding [I (p)] 
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inf„ 1 >o i u 2 >o A + (ui, v 2 )/v\. The argument consists of 3 steps. First, we ar¬ 
gue that inf t , 1>0j j ;2 >o A + (tq ,v 2 )/v\ =mm vl>0 , V2 >oA + (v 1 ,v 2 )/vi. Second, we 
argue that every local minimum is a KKT (Karush-Kuhn-Tucker) point, 
which will be defined shortly. Third, we argue that there is exactly one 
KKT point (ui,U 2 ) G (0,oo) x [0,oo); hence, this KKT point must be the 
global minimum of min. ul> o ii;2 >o h + (v\,v 2 )/vi. 

For v\ > 0 and v 2 > 0, clearly, A + (wi, v 2 ) > —M + (9\,9 2 ) > 0. Hence, A + (tq, v 2 )/v\ 
is bounded below by 0 and explodes as v\ decreases to zero. Furthermore, 
by Proposition 3.1 in [4], A + (ui,U 2 )/fi explodes as the norm of v becomes 
large. Also, for each v\ > 0, we know that h. + (vi,v 2 )/vi is continuous in 
v 2 as v 2 converges to 0. Hence, inf„ 1> o^ 2 >o A + (ui, v 2 )/v\ must be a local 
minimum (with v\ > 0) of 

(4.2) minA + (ui,U2)/ui 


(4.3) 


s.t. v 2 > 0. 


It is known that every locally optimal solution to a constrained, differen¬ 
tiable, nonlinear optimization problem with linear constraints must be a 
KKT point; see 14.37 of [13] or Section 3.5 of [3]. Now, we argue that our 
constrained, nonlinear optimization problem is differentiable. 

The compound Poisson distribution has an infinite support over the inte¬ 
gers; hence, by Proposition 3.1 in [4], there are unique values 9+(v ), 9^(v ), 
such that 


A + (u) = 9f(v)v i + 9^{v)v 2 - M + (9 + {v)). 


(4.4) 


To explore the relationship between v and 9 + (v), fix v and find 9 + that 
maximizes Q\V\ + 9 2 v 2 — M + (9i,9 2 ). Taking derivatives and setting them 
equal to zero yields 



Thus, 9f(vi,v 2 ) and 9^(vi,v 2 ) determine v\ and v 2 . Furthermore, this map¬ 
ping from 9 + to v is a smooth bijection. Since the Jacobian, which is positive 
definite since it is also the Hessian of the strictly convex M + , has a nonzero 
determinant, it follows from the inversion theorem (see [2]), that 9+(v\,v 2 ) 
and 02~(ui,U2) are smooth functions of {v\,v 2 ). Hence, our nonlinear pro¬ 
gramming problem is differentiable. 

For our problem, (v\,v 2 ) is a KKT point if there exists a Lagrange mul¬ 
tiplier u with 


(4.7) uv 2 = 0 


(complementary slackness), 
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(4.8) «>0 


(sign restriction), 


(4.9) 


(u) = V( a+ (wi,u 2 )M) 


(gradient equation), 


(4.10) u 2 > 0 


(constraint). 


Among other things, the proof of our next result shows that any KKT point 
(v\,v 2 ) £ (0,oo) x [0,oo) must have u > 0; hence, complementary slackness 
implies that v 2 = 0. Consequently, there is no “best” path in Fp. that is, the 
inf pe F i I(p) is not attained in Fp 


Theorem 1. For the paths in the interior, we have 
inf [/(p)] = inf A+ (ui, v 2 )/v 1 = 0 \, 

p&Fi > 0,112 >0 


where 9 b = 9 + (v b ) and u 6 = (uj,^) is the unique KKT point in (0, oo) x 
[0, oo) point for (4.2) and (4.3). The point u* is the unique solution to 


(4.11) uj> 0, 

(4.12) v b = 0, 

(4.13) M + (9 + (v b )) = 0. 


Equivalently, but more usefully, 6 b is the unique solution to 


(4.14) 

(4.15) 

(4.16) 


0 \> 0, 

d[M+{e\,e b 2 )\ n 
d0 2 

M + {6 b ) = 0, 


Proof. First, 


(4.17) 


<9[A + (vi,v 2 )\ 
dvi 


0i"(vi,v 2 ), 


which can be seen by starting with (4.4) and using the left most equation 
in (4.5) and (4.6). It follows that 


(4.18) ^ ~ = l v ^t( v D^ 2 ) “ A + (ui ,v 2 )\/v\ 

and 

(4.19) 9[AR^^)M] = + 

ov 2 

In the remainder of this proof assume that (vi,v 2 ) is a KKT point, and 
we know that there exists at least one such point. Thus, using (4.19), the 
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Lagrange multiplier is u = Q^iyi, v 2 ), and using complementary slackness, 
(4.18) becomes 

a[A>i ^)M =M+{e+(ir) )/ vl 

OV 1 

Now, we argue that u > 0. Assume the contrary, that is, that u = 0. 
Thus, the gradient of A. + (v\,v 2 )/v\ at (vi,v 2 ) is zero, which means that 
M + (6 + (v)) =0 and (? 2 ~(ui, u 2 ) = 0. Solving M + ( 6 f, 0) = 0 yields exp(6 , J t ") = 
Mi/(Ai +At 2 r 2 ,i). 

Substituting into (4.5) and (4.6) yields 

vi(0f , 0) = fi\ (Ai + /x 2 r 2> 1 ) 


and 


fa + t , (Ai + /i 2 r 2 1 ) 

v 2 Wi ,0) = A 2 - ir 2 r 2j0 + Mm ,2 - 

Mi 

= A 2 - /i 2 r 2;0 + ri )2 Ai + M 2 ?" 2 ,m ,2 - /r 2 r 2j i 


9i 

92T2,1— j-- 

Ml (Al+M2L2,l) 

91 


(Ai +// 2 r 2j i 


= (A 2 + Al7T j2 ) - M2(l - ^ 2 , 1 ) + M2^2,l r !,2 - 92l"2,l- 


91 


(Ai +/i 2 r 2 ,i) 


(Al + M2^2,l) 

= A 2 (l - ri >2 r 2) i) - M 2 (l - n, 2 ^ 2 , 1 ) + M 2 ?" 2 ,i ^1 - 
= (A 2 - M 2 )(1 - Tl,2^2,l) + 92 T 2,1 ^1 - 
= (A 2 - M 2 )(l - r lt 2 r 2 ,i) - 92 r 2 ,i 


Mi 


+ M2^2,l) 

vi(9f ) 0) 

(Ai T 9 2 r 2 ,i) 


< 0 since A 2 < /i 2 , v\ > 0, 


which is a contradiction. Hence, u > 0, and from complementary slackness, 
we have v 2 = 0. 

Equations (4.12) and (4.13) follow from the gradient equation, and are 
rephrased in terms of 6 h in (4.15) and (4.16). Equation (4.14) follows from 
(4.11) since (4.18) must equal 0. Thus, the two sets of equations are equiva¬ 
lent. To show uniqueness, consider the second set of equations and recall that 
M + is a strictly convex function with M + (0,0) = 0. The set of all (0i,0 2 ) 
such that M + ( 01 , 62 ) = 0 is the boundary of a strictly convex set, the egg 
shaped region in Figure 2, containing (0,0) and (0i,0 2 ) > (0,0). There are 
exactly two points on the boundary of this convex set that are tangent to 
vertical lines, that is, satisfy (4.15), but only one of the two, the eastern 
most point on the boundary, satisfies (4.14). □ 
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4.2. The jitter path on the x-axis In this section we consider paths that 
bounce along the x-axis; that is, we consider paths in F that initally have 
the form p(s) = (uis,0) until hitting (1,0) at time T = l/v\, where v\ > 0 
and after time T, the path follows the natural drift path until reaching the 
origin. Let F x C F denote the set of all such paths, which will be the jitter 
paths and the bridge path. To analyze these paths, we view the x-axis as a 
flat boundary as in Definition 8.7 in [15] and IT as a flat boundary process. 

Define the log moment generating functions of the compound Poisson 
process associated with jumps on the x-axis, 

M- (9 1 , e 2 ) = Ai(e 01 - 1) + A 2 (e 02 - 1) + ^n,o ( e - 01 - 1) + ^n, 2 (e 02 " 01 - 1). 
The associated local rate function is 


A (vi,v 2 ) = sup( 6 * 1 X 1 + 02^2 — M ( 61 , 62 )). 
01,02 


Using the same arguments as in the previous section, there exists a unique 
pair 6 ~(v) such that 


A (vi,v 2 ) = 6 X (vi,v 2 )vi+6 2 (v 1 ,v 2 )v 2 -M (6 (v)). 

The local rate function for the path with velocity v = (ui,0) is given by 
A*(ui,0)= inf (/3A + (u f ) + (1 - /3)A~(v~)), 

0<P<l,/3v++(l-p)v^(vi,0) 


which is a good rate function; see (v) of Lemma 8.20 of [15]. Intuitively, 
the path is a mixture with (3 representing the proportion of time above the 
x-axis, while 1 — /3 is the proportion of time on the x-axis. The bridge path 
has (3 = 1; jitter paths have /3 < 1. Finally, we must calculate infpg^ [I(p)]; 
that is, 


(4.20) 


inf 

vi >0,0</3<l,/3v++(l—(3)v~ =(vi ,0) 


f(/3,v 1 ,vf,v^,v 1 ,v 2 ), 


where f(P,vi,vf ,v 2 ,v 1 ,v 2 ) = (/3A + (v + ) + (1 — (3 )A (v ))/v\. Our argu¬ 
ment will be similar to the last section. First, we argue that (4.20) equals 


(4.21) 


min 

vi>0,0<(3<l,(3v++(l-/3)v- 


, J(P,vi,vf,v 2 ,v 1 ,v 2 ). 

=(V1,0) 


Just as for A + , A~(v~) goes to infinity as |u - [ diverges. Furthermore, 
A~(v~) goes to infinity as v 2 | 0 and is infinite if v 2 < 0. Hence, / goes 
to infinity as f3 [ 0. If [3 > 0, / goes to infinity as vi j 0. Hence, (4.20) must 
be a local minimum of 


(4.22) minf(P,v 1 ,vf,v^,vj;,v 2 ) 

(4.23) s.t. g l (P,vi,vf,v^,vj;,v 2 ) = f3<l, 

(4.24) g 2 ((3,v 1 ,vf,v£,vi,v 2 ) = -v 1 + (3v^ + (1 - f3)vjj = 0, 

(4.25) g 3 (P,v 1 ,vf,v^,vjj,v 2 ) = f3v^ + (1 - f3)v 2 = 0. 
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To show that this constrained, nonlinear optimization problem is differen¬ 
tiable is almost identical to the argument in the previous section includ¬ 
ing the argument that 9~(y 1 ,^ 2 ) is smooth. The constraints are no longer 
linear since there are terms like (3vf. However, the gradients of the three 
constraint equations are linearly independent; hence, from 14.37 of [13] or 
Section 3.5 of [3], all points are regular and every local minimum must be a 
KKT point. The remainder of the argument is to determine the KKT points 
(/ 3 , v\, vf, V 2 , Vi , V 2 ) £ (0,1] x (0,oo) x (—00,00) x (— 00,0 ] x (—00,00) x 
[0, 00). Note that we did not include > 0 and < 0; these will be sub- 
optimal since A - (uf , ) = 00 when < 0. 

For our problem, ((3, v\, vf , v£, iff”, i^") is a KKT point if there exists 
Lagrange multipliers ui, U2 and 113 satisfying 


=0 

u\ < 0 
3 

^2u i Vg i = X / f 
1=1 

v\ + /3vf + (1 - (3)vi = 0 
P v 2 + (1 - P)V2 = 0 


(complementary slackness) 
(sign restriction), 

(gradient equation), 

(constraint 1), 

(constraint 2), 

(constraint 3). 


Now assume that ((3,vi,v+,v£,v± ,v%) is a KKT point. The bottom four 
components of the gradient equation imply that v\U 2 = df(v + ) = Q^(y~) 
and V 1 U 3 = ^ 2 ^(u + ) = Q^iy-). Using these in the second component of the 
gradient equation implies that [3M + (9 + (y + )) + (1 — P)M~(9 + (v + )) = 0. 
Since v\ > 0, the second component of the gradient equation also implies 
that U 2 > 0, which means that v\U 2 = ^(u"*") > 0. The first component of 
the gradient equation reduces to u\ = — [M + (9 + (v + )) — M~(9 + (v + ))]/v±. 

First, consider the case when j3 < 1. By complementary slackness, u± = 0 
implying that 


(4.26) M + (0 + (F f )) =M _ (6 » + (F f )) = 0. 


These level sets enclose convex regions which intersect at (0,0) and possibly 
other points (#i~,# 2 ~)- However, since v£ = <>A1 ^ ~ must be (strictly) nega¬ 
tive, we are only interested in solutions to (4.26) in the lower portion of the 
egg in Figure 2; that is, going clockwise from (but not including) 9 b along 
M + = 0 to the other solution of (4.15), which is the western most point of 
the egg shaped region defined by M + = 0. Since 9f (u + ) > 0, we can further 
restrict the region to the arc going clockwise from the 9 b to the origin. 



16 


R. D. FOLEY AND D. R. MCDONALD 


By calculation the points [besides (0,0)] where M + (0) = 0 and M~(9) = 0 
cut the 0i-axis are (log(/ii/(Ai + P 2 i" 2 ,i)), 0 ) and (log(/x*/Ai), 0 ), respec¬ 
tively. Since the first coordinate of the latter is positive and greater than 
the first coordinate of the former, it follows that 9 2 > 0 if 0 ^ > 0 ; thus, 
we can restrict attention to the arc going clockwise from 9 h to the 0i-axis 
at (max(log(//i/(Ai + 1 ^ 2 f 2 , 1 )), 0), 0). There can be at most one such point 
since ^ > 0 . If there is such a point, label this point as 0 J = (9{,9%); 

otherwise, define 9 3 = 9 h . Note that there can be no other KKT points with 
(3 < 1. To see that 0 J determines a KKT point, recall that vf = —^ and 
v\ = ^ for i = 1,2 with vf > 0, v 2 <0 and v 2 > 0. Since (1,0) lies in 

the convex hull of v + and v ~, there exists a unique f3 and v\ with 0 < (3 < 1 
such that (3v + + (1 — (3)v~ = (ui,0). Let u\ = 0, 112 = 9{/v\ and U 3 = 9 3 2 /vi. 
It is straightforward to show that these values satisfy the KKT conditions. 

To complete the first case with (3 < 1, we need to argue that there cannot 
be another KKT point with (3=1. If there were another such KKT point 
with (3 = 1. then our optimization problem is a special case of the one in the 
previous section; hence, the solution must correspond to 9 h , in which case the 
first component of the gradient equation would imply that u\ = M~(9 b )/v 1 . 
However, this would violate the sign restriction on u\ since M~(9 b ) would 
be strictly positive when 0 J 7 ^ 9 b . Consequently, when 0 J 7 ^ 9 b , we have a 
unique KKT point, which must be the global minimum. 

Suppose there is no solution 0 + to (4.26) on the clockwise arc following 
M + = 0 from 9 b to the 0i -axis. Then there cannot be a KKT point with (3 < 
1. However, we know that there is a global minimum, which must be a KKT 
point. Hence, any KKT points must have (3 = 1. If (3 = 1, the optimization 
problem in this section becomes a special case of the optimization problem 
in the previous section, and it follows that there is a unique KKT point 
corresponding to 9 b . 

Theorem 2 . Among the jitter and bridge paths, we have 
( 4 . 27 ) inf [I(p)} = 9{. 

p£F x 

Furthermore, the following conditions are equivalent: 

1. 9 3 <9 b , 

2 . a jitter path is optimal in F x , 

3. (3 <1. 

4. M~(9 b )> 0 and 

5. p <1, where 
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Proof. To show (4.27), note that 9 3 determines the global optimum, so 
we need only evaluate the objective function at the point determined by 9°. 
If = 0 b , the objective function becomes identical to that of the previous 
section, which we know evaluates to 9\=9\. If 0 ] < 6 b , the objective function 
simplifies to 9 

The first three conditions are clearly equivalent from the discussion prior 
to the statement of the theorem. For the fourth condition, it suffices to show 
that the directional derivative of M~ at the origin tangent to the level curve 
M + = 0 and going counterclockwise is strictly negative; that is, M~ = 0 
pierces the egg M + = 0. Since M - (0,0) = 0, it would follow from continuity 
that the level curves for M + = 0 and M~ = 0 would have to intersect along 
the segment M + = 0 somewhere between 9 b and the origin going clockwise 
from 9 b . (From the proof of the previous theorem, we could further restrict 
attention to the segment lying in the first quadrant.) 

If ( x + ,y + ) denotes the gradient of M + at the origin, then (— y + ,x + ) is 
tangent to M + = 0 at the origin and points in a counterclockwise direction. 
If ( x~ ,y ~) denotes the gradient of M~ at the origin, then the directional 
derivative of M~ at the origin tangent to the level curve M + = 0 and go¬ 
ing counterclockwise is — y + x~ + x + y~ = // 2 Ai — y\\irip + H \^2 — P 1 P 2 ~ 
Hi A 2 + H 2 ^ 2 i" 2 ,i + f 1 ,+ i"i, 21 ^ 11 ^ 21 ^ 2 , 1 - This directional derivative is a 
strictly decreasing function of Hi since the derivative with respect to Hi is 
H 2 ( 1 - n,2^2,i)(/92 - 1 ) < 0. If Hi = (Ai P 2 ^i) /(1 P 2 ); then the directional 
derivative is zero and M + = 0 and M~ = 0 are tangent at the origin. How¬ 
ever, from Proposition 1, stability requires that Hi > (Ai — P2^i)/(1 ~ P 2 ), 
which ensures that is decreasing in the direction (— y + ,x + ) and pierces 
the egg M + = 0. 

To complete the proof, using (4.6), it is straightforward to show that the 
fifth condition is equivalent to uj (9 J ) < 0; recall that (9 b ) = 0 and see Fig¬ 
ure 2. 

□ 

Corollary 1. For the Jackson network with h\ = Hi> a jitter path is 
better than the bridge if (2.2) holds. 

Proof. In the Jackson case exp (0{) = pf l and exp^) = {r 2 ,o+i' 2 ,ipf 1 ). 
Using this, p < 1 by substitution and Theorem 2 gives the result. □ 

Now, we will try to derive an explicit expression for 9 ] by locating the 
points 9 / (0,0), where M + (9) = M~(9 ) = 0. Solving M~(9) = 0 yields 

Ai + A 2 + Hi ~ ~ P*i r ifl e ~ 61 

A 2 + Pi r i,2&~ ei 


(4.28) 


exp ( 02 ) 
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After substituting this into M + (0) = 0 and simplifying, we see that x = 
exp(9j) must be a positive solution to the quadratic equation 

(4.29) ax 2 + bx + c = 0, 

where 

a = — /d)(A 2 + Ain^Ai — A2/^2(A2?"2 ,i + Ai), 

b = —(n i — m)(A2 Mi + MiAri,2 + A 2 + Afri t 2 + 2AiA2r%2 + Ai^rqo) 

+ /^A 2 At2(l - 2r-i )2 r- 2 ,i) - MiAiM2n,2, 

c=(nl~ /ri)^ri i0 (A 2 + Airi j2 ) + (mi) 2 M2?T,2^1,0 + 2^2,0)- 

Given such a solution x, then 6 \ = log(x) and 62 is determined from (4.28). 
Now that we have a point 0 where the level curves intersect, we need to 
determine if it lies on the arc M + = 0 going clockwise from 6 b to the 01- 
axis. Since 0i > 0, we are only interested in x > 1. If, in addition, uj (0) < 0, 
then 0 = 6 : . Note that there can be at most one x giving a solution that 
satisfies these conditions. If there is no such solution x , then 0 jf = 6 b . 


Remark. The paper [10] provides an alternative to the Schwartz-Weiss 
flat boundary approach to determining the local rate function for the jitter 
and bridge paths. Rather than mixing the two vectors and v~ with the 
weights P and 1 — /3, [10] simply expresses A*(u,0) = sup 7 (7 • v — A(y)), 
where A(0) = log(r(J 7 )), r(J 7 ) is the spectral radius of J 7 and J 7 is the 
Feynmann-Kac transform of the kernel J of the Markov additive process 
associated with the modified Jackson network as defined in [9]. If we could 
show that the mapping from 7 to v is a smooth bijection, then since A( 7 ) 
is convex, we could represent 


inf 

w>0 


A*(u,0) 

v 



Taking derivatives, the minimum occurs when 


A( 7 ) 


d 2 A( 7 ) / 

dj 2 / 



= 0 . 


Since A( 7 ) is convex by Lemma 2 in [10] and since c g)) = v > 0, it follows 
that the minimum occurs when A ( 7 ) = 0 and the large deviation rate is the 
associated 7 . This just means that the large deviation rate is the choice of 
7 which sets r(J 7 ) = exp(A( 7 )) = 1, which agrees with the results in [9]. 

The problem is the assumption that the mapping from 7 to v defined by 
sup 7 (7 ■ v — A(y)) is smooth. Lemma 3.3 in [12] shows that if 7 E U r . where 


U r = { 7 : (p( 7 , A) = 1 for some A = A( 7 ) < 00 }, 
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and if ( 7 , A(t)) 6 W, as defined in [12], then A(-) is differentiable at 7 . Unfor¬ 
tunately, we are particularly interested in cases when 7 may not belong to U r 
because these give rise to bridges. It appears that in the two-dimensional 
case differentiability may follow from the explicit description of A ( 7 ) (private 
communication with Ignatiouk-Robert), but in n dimensions this is far from 
clear. Hence, at first blush, it appears that the optimization problem would 
be more difficult to solve than the optimization problem arising from the 
Schwartz-Weiss approach where A* is represented as a convex combination 
of smooth pieces. 

Another advantage of the Schwartz-Weiss approach is that it distinguishes 
between the bridge (/3 = 1) and the jitter (/? < 1) paths. It is not clear 
that the approach in [10] makes this distinction. On the other hand, the 
results in [ 10 ] are much more general since they apply in n dimensions and 
apply to permeable or impermeable boundaries. Note that the condition 
M+{ 6 b ) >M~{ 6 b ) in Proposition 10 in [10] inspired condition 4 in Theorem 
2 . 


4.3. Cascade paths climbing the y-axis. Finally, we consider fluid paths 
that go up the y-axis to a height (0 ,h) and then go down and across ( 1 , 0 ). 
For h > 0, we will refer to these paths as cascade paths since the customers 
build up in the second queue and then “cascade” into the first queue for the 
large deviation. The path with h = 0 is none other than the bridge path. 
Note that a bridge path up the y-axis followed by a cascade is not optimal 
because this would give a nonlinear large deviation path in a domain with 
constant jump rates. Consequently, the least action path jitters up the y-axis 
and then cascades. Let F c C F be the set of cascade paths and the bridge 
path. 

We wish to find conditions for a jitter path along the y-axis when p 2 < 1, 
but pi but may be greater than one. The investigation of jitter paths in 
Section 4.2 can be used provided we interchange the x- and y-axes. Let 
6 = { 61 , 62 ) be the solution analogous to 6 J in Theorem 2. Thus, equivalent 
to (4.26), we have 

M + { 6 iA) = M{ 6 i,~ 6 2 ) = 0, 

where 

M( 6 i, 6 2 ) := Xi ^ 1 - 1) + A 2 (e 0 ' 2 - 1) +p 2 r 2 Ae - §2 - 1) - 1). 

Subtracting M from M + gives exp(<?i) = r^o + r\p exp($ 2 )- Substituting 
this into M{6) = 0 gives exp(# 2 ) equal to 1 or / 12 /A 2 , and the former solution 
can be eliminated since it corresponds to 6 ± = 62 = 0. In order to have a jitter 
path, condition 4 in Theorem 2 must hold; that is, 

Xie Sl + P2 , r2,ie~^ +Sl 

p :=-=-=—=- < 1. 

hi r ifi e ~ dl + piripe 02 - 01 
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Substituting exp($ 2 ) = P 2 /A 2 gives p = \\ exp($i)/pi < 1; that is, (2.3) holds. 

We conclude there is no cascade path unless (2.3) holds, but if it does, 
then the large deviations rate of paths from (0,0) to (0 ,h) is /ilog {pf 1 ), 
which can be seen by simplifying (4.29) when = p\ to obtain 9{ = pf 1 . 
Moreover, the large deviation rate of paths from (0, h) to (1,0) is 


inf 

D>0 


A + (u, —vh) 
v 


After (1,0), the process follows the natural drift path back to zero. Thus, if 
path p e F c reaches height h, we have I(p) = hlog(pf 1 ) + inf„>o A ( v e~ L,h \ 
We wish to find inf pg j7 c /(p) or, equivalently, inf/ l >o jt) >o f{h , v), where f(h, v ) = 

hlog(pf 1 ) + A ( v ’~ vh ^ . Since / is continuous, positive and diverges as v j 0 
or as v —> 00 or as h j 0, we know inf/ l >o ! t,>o f(h, v ) must be a local minimum 
of 


(4.30) min f(h,v) 

(4.31) s.t. h > 0. 

This constrained nonlinear optimization problem is differentiable with a 
linear constraint so every local minimum must be a KKT point. In order for 
(h, v ) to be a KKT point for this problem, there needs to be a corresponding 
Lagrange multiplier u so that h,v and u satisfy the following: 

(4.32) uh = 0 (complementary slackness), 

(4.33) u>0 (sign restriction), 

( 4 - 34 ) ( 0 ) = V/(/i,v) = ^ ) (gradient equation), 

(4.35) h > 0 (constraint). 

As in the previous sections, there will be a unique KKT point (h,v) E 
[0, 00 ) x (0, 00 ). 


Theorem 3. Define 9 C = 9 + (v, —vh), where (v, h) minimizes (4.30) sub¬ 
ject to (4.31). It follows that 

(436) inf [I(p)} = 9 c v 

p&f c 

If pf 1 > + ri 5 2 P 2 1 an d < ® 2 > then the minimum action path in 

F c is a cascade of height h > 0 with 9 C = (log(pj" 1 ),log(p^" 1 )); otherwise, the 
minimum action path in F c is a bridge with 9 C = 9 b . 


Proof. We will argue that 9 C corresponds to the unique KKT point 
in (h,v) E [0, 00 ) x (0,oo); hence, this point must correspond to the global 
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minimum. Note that the second component of the gradient equation implies 
that we must be looking for a solution corresponding to a point 9 + with 
M + {9 + ) = 0. This point 9 + determines v, h and u since v = vf(9 + ), h = 
( 9 + )/v, and u = log^ 1 ) - 6 %. 

In order to have a KKT point 9 + with u > 0 we would have to have 
h = 0, which corresponds to a bridge path. Hence, 9 + = 9 b and by (4.34), 
this occurs if and only if (2.3) is false or log(p^ 1 ) > 0 2 - (4.36) follows by 
substitution. 

Alternatively, in order to have a solution with u = 0, we would have 
to have h > 0, which corresponds to a cascade path. Hence, 9 + = 9 C and 
this occurs if and only if (2.3) is true and log(p^ 1 ) < 0 2 . Next, (4.34) re¬ 
quires 9 2 = log^ 1 ) < 0 2 and M + (9 + (v, —vh)) = 0. Moreover, VM + (0 + ) = 
(v,—vh). The two solutions to M + (9 + (v,—vh)) = 0 and 9£ = log^^ 1 ) are 

(logCpr 1 ),!^^ 1 )) and (l°g( r ll0 + r ll2 P2 1 ) do S(/ 3 2 _1 ))- 
Substituting the second solution into (4.5) and (4.6) gives 

VM + (0 + ) = ((Ai + A 2 r 2 ,i)e 0 i + - pie~ 0 +(ri jO + p 2 Vi j2 ), 

A 2 p 2 1 — A 2 r 2j o + pi?T i2 p 2 1 e dl — A 2 r 2j ie e ! ). 

Notice that the first coordinate equals Aie^i" — pi = Ai((ri j o Tripp^" 1 ) — pj" 1 ). 
This is negative if (2.3) is true and this is impossible since VM + (0 + ) = 
(v, —vh) with v > 0. 

Substituting the first solution into (4.5) and (4.6) gives 

VM + (0 + ) = (pi - Airpo - Am )2 p 2 :1 ,P 2 - A 2 r 2)0 - A 2 r 2j ipf 1 ) = (v, -vh). 

The first component is Ai(pj) 1 — (rqo + n,2P 2 1 )) and this is positive since 
(2.3) is true. The second component of VM + (0 + ) is negative since 9 + is on 
the level curve M + {9) = 0 and 9% < 0 2 . (2.2) fails. Equation (4.36) follows by 
substitution. 

□ 

4.4. Summary of flat boundary approach. Combining yields the following 
result. 

Theorem 4. For the modified Jackson network, 
lim l log P(W € F(£)\W(0) = (0,0)) 

a —KX) t 

= lim ilog P(W t G F(l)\W e (0) = (0,0)) 

£—>00 Z 
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If9\ < min{log(p 2 1 )? ^ 2 }? then minimum action is 9\ and the minimal action 
path is a jitter path along the x-axis. If log(p 2 x ) < min{# 2 ,0 2 }, then the 
minimal action is Of and the minimal action path is a cascade path that 
initially climbs the y-axis. Otherwise, 0\ = 9{ = 9(, the minimal action is 
9\, and the minimal action path is a bridge. 


Proof. The only part of the proof that is not straightforward is showing 
that if log (p^ -1 ) < min{# 2 ,0 2 }, th en the condition pf 1 > r 10 + r\gp 2 l of 
Theorem 3 holds automatically. To prove this, we will show that if pf 1 < 
r i,o + r ipP 2 1 and log(p 2 X ) < $ 2 > then 9 3 2 < log^pf 1 ). In other words, a jitter 
path will be the minimal action path—not a cascade. 

The (convex) function g~ defined by M~ {9\, g~ (# 1 )) = 0 is given at (4.28). 
Simplifying, 


g (6>i) = l + 


pt(l-e~ ei ) + Ai(l-e ei ) 
A 2 + g\ri,2e~ ei 


Differentiating g~ with respect to g\ gives (1 — exp(— 0 i))/(A 2 + Pi^i ,2 x 
exp(—#i )) 2 and this is strictly positive for 9\ G [0,log(/i*/Ai)], where the 
endpoints are the zeros of g~. Thus, g~ is strictly increasing in g* on 
[0,log(//*/Ai)]. Furthermore, lim /i ._ voo g _ (log(ri i o + J'i i 2 P 2 ' 1 )) = log^ 1 ), which 
is a point on M + = 0. If pf 1 < rig + r\ )2 pf 1 , then the easternmost solution 
9 to M + (9) = 0 with 9 2 = log^ 1 ) is (log(ri j0 + rigpf 1 ),\og{pf 1 )). It fol¬ 
lows that M + = 0 and M~ = 0 intersect at a point in the positive quadrant 
going clockwise from (log^^o + n, 2 P 2 1 ))l°g(P^ 1 )) along M + = 0 before hit¬ 
ting the axis. If, in addition, log(p^" 1 ) < # 2 , then the point of intersection 
defines (9 J , which must be the minimal action path. 

There is one boundary case which also must be eliminated. We must 
show it is impossible that 9 J = 9 C < 9 b . If this did happen, then 9 J = 9 C = 
(log(p^ 1 ),log(p 2 x )). Since lies on M + and M~ , it follows that ^“(log^ 1 )) = 
pf 1 . Since 9 J < 9 b , it follows that p < 1 and this reduces to p 2 1 < r 2 g + 
ra.ipr 1 - Now notice that when g\= g\, that is, in the Jackson case, we 
have ^“(log^j) 1 )) = r 2 g + r 2 gpf l and this is strictly greater than pf 1 as 
we just showed. This eliminates this case if g\ = g\. Moreover, we showed 
above that g~(9\) is increasing in g*, so it will never be possible to solve 

p-(io g (pr 1 ))=P2 1 . □ 


Corollary 2. For the Jackson network with g\ = g\, if (2.2) fails and 
(r 2t o + r 2 gpf 1 ) > pf 1 , then a cascade path is optimal, whereas if (2.2) fails 
and (r 2 g + r 2 g P\ l ) = P 2 l , then a bridge path is optimal. In both cases the 
minimal action is log(pj) 1 ). (This is no surprise since the steady state ir is 
of product form even in the cascade or bridge cases.) 
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Proof. In the Jackson case exp(ffj) = p 1 1 and exp (# 3 ) = (^ 2 ,o+^ 2 ,iPi X )- 
Hence, if ( 72,0 + t^i/rf 1 ) > p^ 1 , then exp(# 2 ) > P 2 1 ■ Using this, p > 1 by 
substitution and Theorem 2 shows = 6 b . This means exp (# 3 ) > P 2 } and, 
moreover, if (2.2) fails, then (2.3) must hold. By Theorem 3 it follows the 
minimum action path is a cascade and Of = log(pj) 1 ). 

If (t 2,0 + f 2 ,iP\ 1 ) = P 2 1 ) then by the above argument exp(# 2 ) = p^ 1 . 
Hence, by Theorem 3, the bridge path is optimal and 6 \ = 0{ = log(p)” 1 ). 
□ 


5. Another model exhibiting the bridge phenomenon. We can use the 

above methods to obtain the rough asymptotics of n(£,y) for a fork network 
as introduced by Flatto and Hahn [7]. This model was later extended and 
described as the bathroom problem in Chapter 16 of [15] . Couples arrive at 
a cinema according to a Poisson process with rate v and immediately visit 
the men’s and ladies’ room. The service rate the men’s queue a, while the 
rate at at the lady’s queue is (3. There are also separate arrival streams, with 
rate p for single women and rate A for single men. 

We are interested in a large deviation of the men’s queue, so let this be 
the first queue and the ladies’ queue the second. We could use the above 
methods to obtain the rough asymptotics of n(£,y), but we don’t have to 
because the work is already done in Section 16.2 in [15]. In the (unrealistic) 
case p + av/(\ + v) < (3, they show the most likely path for the men’s queue 
to reach a high level £ before returning to zero is a jitter path along the x-axis 
and the large deviation rate is a/(\ + v) in agreement with [11]. In the (more 
realistic) case when p + a.v /(A + u) > f3, the results in [15] show the most 
likely path for the men’s queue to reach a high level l before returning to 
zero is a path through the interior (again in agreement with [11]). Moreover, 
among jitter paths, the large deviation rate (for the more realistic case) is 
minimized by paths spending zero time on the boundary. 

This just means that if, in addition to the men’s queue getting big, we 
require that the ladies’ queue remain small, then the large deviation path is 
a bridge path. The exact asymptotics of 7 r(£, y ) are calculated in [9] and are 
found to agree with those in [7]. The derivation of these sharp asymptotics 
confirms our intuition about the bridge behavior. 
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